Exploring Exotic Superfluidity of Polarized Ultracold Fermions in Optical Lattices 
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Enormous interest has been paid to ultracold 
Fermi gases due to the interplay between Cooper 
pairing and strong correlations [11121, Isl]. Beautiful 
experiments on the superfluidity have been per- 
formed in these systems with unequal spin popu- 
lations [4, y, Isl, III, Isj, IsJ . Arrest ingly, it was found 
that the superfluid paired core is surrounded by a 
shell of normal unpaired fermions while the den- 
sity distribution of the difference of the two com- 
ponents becomes bimodal. Here we explore theo- 
retically the novel superfluidity of harmonically- 
trapped polarized ultracold fermionic atoms in a 
two-dimensional (2D) optical lattice by solving 
the Bogoliubov-de Gennes equations. The pair- 
ing amplitude is found to oscillate along the ra- 
dial direction at low particle density and along 
the angular direction at high density. The for- 
mer is consistent with the existing experiments 
and the latter is a newly predicted Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state, which can be 
tested in experiments. 

Experimentally, optical lattices artificially created by 
standing- wave laser fields [l^, [Hi 12 1 provide unprece- 
dented experimental tunability, such as a precise control 
of model Hamiltonian parameters including the interac- 
tion strength, lattice geometry and number of atoms. 
This versatile tool opens a door to explore rich and novel 
physical phenomena in various strongly correlated sys- 
tems. Experimental exploration of the novel superfluid 
phase of imbalanced fermions in optical lattices is ex- 
pected to be conducted in the near future, while in the 
theoretical aspect, there exists a long-standing debate on 
the nature of the exotic pairing state for the imbalanced 
systems. The inhomogeneous superconducting state, 
known as the FFLO state, was predicted about forty 
years ago, which contains nonzero momentum Cooper 
pairs with the superconducting gap exhibiting a periodic 
real-space modulations [3] ■ There have been suggestions 
for the possible FFLO phase for some of the high- Tc, lay- 
ered organic, and heavy fermion superconductors at very 
low temperatures and high magnetic fields [3, • This 
subject has also been studied in nuclear and high-energy 
physics lldl . Despite of many previous literatures in this 
field ^[M, E M, [ii, Hi, Q 0, S Hi Hi , it is still 



awaited to explore theoretically and unambiguously the 
exotic FFLO state with imbalanced ultracold fermions in 
optical lattices. 

In this letter, we present an unrestricted Hartree-Fock 
study of a two-component fermionic atoms in a harmoni- 
cally confined 2D optical lattice. Due to the competition 
between the ferromagnetic and s-wave superfluid phases 
in the presence of harmonically trapped potential, the 
compromised ground state may exhibit rich patterns of 
order parameter configurations. Our main findings are 
as follows, (i) In the low fermionic filling regime, the 
superfluid-pairing gap modulates its sign along the radial 
direction, indicating an FFLO state. The calculated den- 
sity profiles of the two components of the fermions show 
a striking resemblance to the experimental observations, 
(ii) In the high filling regime, a superfiuid-pairing ring 
appears. As the imbalance increases, the pairing gap 
oscillates along the angular direction, signaling a novel 
FFLO state, (iii) A variety of intriguing distributions of 
the order parameter, such as square lattice FFLO state, 
emerges in the intermediate filling regime. In addition, 
we also observe that, as the imbalance increases, the su- 
perfluidity is suppressed and vanishes at a critical imbal- 
ance amplitude, implying the existence of the Clogston 
limit [ig. 

To simulate the neutral fermionic atoms in a 2D optical 
lattice, we begin with an effective tight-binding Hamil- 
tonian, which captures the basic interplay between fer- 
romagnetism and s-wave fermion-pairing in the presence 
of a confining potential, 

H = -tY^ c^^Cj,, +^(ei (1) 
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where Ci„ is an annihilation operator for an atom at site 
i (at position r,) with spill fj and Ti^fj — ^ifj^ia is 

the 

number operator, t is the effective nearest-neighboring 
hopping integral. The spin-dependent chemical potential 
due to the Zeeman splitting is defined as /lo- /i — ah. 
Ci = ^muP\Y\ — rop is the harmonic confining potential 
at site i, and ro is the position of the trap center. 

At the mean-field level, the model Hamiltonian can 
be diagonalized by solving the resulting Bogoliubov-de 
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Gennes (BdG) equations self-consistently S^, 31 1, 
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where i/ij,o- 
Hamiltonian, A 



- (e; + ^J■a■)Sii is the single particle 
— Ai(Sii is the pairing order parame- 
ter, and (uf^jVfg) are the Bogoliubov quasiparticle am- 
plitudes at the site i. The atom density and the s-wave 
pairing order parameter satisfy the self-consistency con- 



ditions, 



and A, = l^E„«T<i* + <T<r)tanh(i;„/2/csr), re- 
spectively, where T is the temperature. The local magne- 
tization is defined as nii = riii^ — , which can be tuned 
by changing the magnetic field in the present model. In 
a laboratory optical lattice, can be tuned directly al- 
though there is no real magnetic field. In our calculations 
below, we set t = 1 as the energy unit and lattice spac- 
ing a = 1 as the length unit. We consider a very low 
temperature T = 0.001, and the amplitude of confining 
potential \mu}^a^ jt ~ (0.016 - 0.025). The dimension 
of lattice is 42 x 42 with open boundary conditions and 
the trap center is located at (21,21). All the order pa- 
rameters or mean fields are calculated self-consistently 
by solving the BdG equations iteratively, with their ini- 
tial values randomly distributed on the lattice. In the 
case of multi-solutions, we compare their corresponding 
free energies to determine the most energetically favored 
state. 

The applied magnetic field may cause the density im- 
balance of the two species of fermions (spin-up and spin- 
down). The induced ferromagnetism competes with the 
superfluid pairing state, resulting in the emergence of ex- 
otic FFLO phases. In the special case = and in 2D, 
the pairing gap of the ground state of Hamiltonian (1) os- 
cillates in space with a stripe-like pattern for the s-wave 
and a square lattice pattern for the d-wave [s^l- In the 
crossing region of the positive and negative pairing ampli- 
tudes, there exist nodal lines where the ferromagnetism 
appears. The presence of the trapped potential may add 
up to the new complication. In the widely used local 
density approximation (LDA), the spatial dependence of 
observable must follow contours of constant trapping po- 
tential in the same way that a spatially uniform system 
does. Our approach provides a consistent way to explore 
novel states beyond LDA. In particular, our results show 
that the compromising ground state depends crucially on 
the Uc (fermion density at the trap center) . The features 
presented below can be classified into three regimes ac- 
cordingly. 

Low-filling regime (ric < Ij. In Fig. 1, we plot typi- 
cally three spatial profiles of the pairing order parameter 
Ai together with local magnetization rrii as functions of 
magnetic field h. At /i = 0, there is no ferromagnetic 
order and the trapping potential may result in the ac- 
cumulation of electron density around the center. As 




FIG. 1: (Color online) Spatial profiles of pairing order pa- 
rameter Ai [(a), (c), and (e)], and local magnetization rrii 
[(b), (d), and (f)] of Eqn. (1) at a low particle filling. 

= 200 fermions on a lattice of 42 x 42. The parameter 
is ^mijj'^o? /t = 0.016. The magnetic field is (top to bottom): 
^ = 0,0.5,1.5, respectively. 



shown in Fig. l(a)-l(b), A^ reaches its maximum value 
at the trap center and decreases to zero far away from the 
center. These results are consistent with those obtained 
by LDA, where no competing orders show up. As h in- 
creases, the ferromagnetic order may emerge and frus- 
trate the conventional pairing state. In Fig. 1(c), our 
calculations show the appearance of superfluid pairing 
gap modulations along the radial direction, indicating 
an FFLO state. At the edge of the modulation, a con- 
tinuous sign change of the pairing gap corresponds to 
a nodal line (A^ = 0). Meanwhile the magnetic order 
parameter rrii is most remarkable around the ring-like 
nodal line so that its density profile exhibits a bimodal 
structure, as shown in Fig. 1(d). In other words, the mi- 
nority fermions are squeezed into the inner core due to 
the pair condensation in the superfluid phase, while the 
unpaired majority fermions are repelled to the outside 
of the core. When the imbalance increases, the minor- 
ity fermions contribute more effectively to the pairing 
states. This result is in good agreement with the experi- 
ments, and also with the previous theoretical studies on 

H, [1^. 



the FFLO state in a continuous model 18, |19|. As h 
progressively increases, the pairing order may further be 
suppressed and the weak modulation of the pairing gap 
appears away from the trap center while the ferromag- 
netic order shows a plateau-like feature nearby the trap 
center. Furthermore, the superfluidity vanishes at a crit- 
ical magnetic field {he ~ 1-6), which corresponds to the 
Clogston limit. Note that the existing experimental re- 
alization (3, H, S 0, B may correspond to the present 
low-filling case. The density profiles obtained in our the- 
ory exhibit remarkable resemblance to the experimental 
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FIG. 2: (Color online) Density profiles along the diagonal of 
the lattice. Ar is the distance from the trap center. Red line 
denotes the magnetization while the black and green lines 
represent the spin-up and spin-down densities, respectively. 
All the parameters are the same as in Fig. 1 



measurements. 

To show the bimodal structure more transparently, 
we display the density profiles along the lattice diago- 
nal through the trap center. The density distribution 
has a four-fold symmetry, consistent with the underly- 
ing lattice structure. As shown in Fig. 2, the bimodal 
distribution of magnetization clearly shows up at finite 
h. There are three regions: a superfluid core with equal 
densities, a partially polarized shell, and a fully polarized 
region. By increasing h, these bimodal structures evolve 
to have more pronounced amplitude while the separation 
between peaks becomes narrower. As h exceeds a critical 
value, the equal density core disappears and the conden- 
sation fraction also vanishes. In such a case, the ferro- 
magnetic order reaches its maximum at the trap center 
and its density profile displays a plateau-like structure. 

High-filling regime (uc close to 2). In this regime, the 
band is almost fully occupied, corresponding to an insu- 
lating state in a conventional solid state system, around 
the trap center. The pairing order can not survive far way 
from the center because of the low fermion density. In the 
intermediate distance from the center, one may naturally 
expect the formation of ring like fermion pairing state at 
the intermediate range. In Fig. 3(a)-3(b), the profile of 
Ai does show such a topological structure at /i = 0. Upon 
increasing h, the ferromagnetic order may emerge and re- 
sult in the frustration of the competing pairing state. It 
is well known that the lower angular momentum state 
has always the lower energy in the infinite system. How- 
ever, in the present case, due to the highly nontrivial 
interplays among trapping potential, magnetization and 
pairing correlations, it is possible that the exotic pairing 
state with higher angular momentum state becomes ener- 
getically favorable. As depicted in Fig. 3(c), the pairing 
gap oscillates along the angular (ring) direction with al- 
ternating positive and negative signs. In Fig. 3(d), it 
shows clearly four regions with abrupt sign change where 
magnetization shows up. This special configuration low- 




FIG. 3: (Color online) Spatial variations of Ai [(a), (c), (e) 
and (g)], and local ferromagnetic order irn [(b), (d), (f) and 
(h)] at a high filling: — 1200 fermions. Parameters are the 
same as in Fig. 1. 



ers the free energy and allows oscillating fermion pairing 
state to remain stable. By further increasing h, there 
will be more oscillations along the angular direction as 
well as more appreciable magnetization, as displayed in 
Fig. 3(e)-3(f). When h = 1.5, the superconducting order 
is almost fully suppressed while the ferromagnetic order 
forms a shell pattern [Fig. 3(g)-3(h)]. 

An intuitive physical understanding of angular FFLO 
state can be given as follows: the main effect of the high 
filling at the trap center is the formation of superconduct- 
ing shell structure as shown in Fig. 3(a). As is known, 
the FFLO is rather robust in ID [33, [s^. If the supercon- 
ducting shell is sufficiently thin like a ID ring, then FFLO 
tends to emerge. That is to say, the spin imbalance will 
lead to the pairing oscillation along the angular direction 
rather than in the radial direction since the former state 
costs less kinetic energy. This novel angular FFLO state 
can be observed in 2D optical lattices where high filling 
regime can be easily achieved. We have checked our re- 
sults for the d-wave pairing interaction, and the similar 
conclusions have also been reached. It is worth to men- 
tion that the simple LDA fails to get such a novel pairing 
state. Here we need to stress the critical role played by 
trapping potential in the formation of the FFLO states 
along radial or angular direction. The appearance of the 
FFLO state minimizes the relevant free energy, which is 
a direct manifestation of the competition between ferro- 
magnetism and superconductivity. The trapping poten- 
tial provides a confined and inhomogeneous background, 
which complicates the competition, but reveals some of 
the novel physics associated with the FFLO states. 
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FIG. 4: (Color online) Spatial profile of (a) pairing order pa- 
rameter, and (b) magnetization rrii at an intermediate filling: 
A'^ = 400 fermions. The parameter is ^mLu'^a? /t = 0.025 and 
other parameters are the same as in Fig. 1. 



Medium-filling regime (rif. G (l,2)j. In this regime, 
rich patterns of order parameter distribution show up. 
In the absence of the magnetic field, the pairing order 
parameter shows a small concaved structure around the 
trap center, as seen in Fig. 4(a). By tuning the magnetic 
field to h = 0.4, the pairing gap modulation along the 
radial direction shows up and the situation is analogous 
to the low filling case [Fig. l(d)-l(f)]. When h = 0.6, the 
pairing order parameter modulates not only along the 
radial direction but also around the trap center, exhibit- 
ing a square lattice pattern as illustrated in Fig. 4(e). 
Meanwhile, its ferromagnetic order exhibits a rich struc- 
ture where there are many local maxima around the trap 
center which correspond to nodal lines in Fig. 4(f). As 
we further increases the magnetic field to ft. = 0.9, the 
square lattice pattern of FFLO state becomes more pro- 
nounced and the associated ferromagnetism exhibits it- 
self as well [Fig. 4(g)-4(h)]. Note that the square lattice 
FFLO state appears in the infinite 2D system with the d- 
wave pairing [34|. When h becomes quite large, the BCS 
pairing order is fully suppressed and the magnetization 
displays a plateau like structure around the trap center, 
as depicted in Fig. 4(i)-4(j). 

So far, there has been no direct evidence of FFLO 
states in recent experiments since the probe of pairing 
gap distribution is very challenging. We now address 
an important issue on how to detect the peculiar real 
space patterns illustrated in Figs. 3 and 4 experimen- 



tally. First we propose a 2D optical lattice experiment 
with the fermionic atoms in the high filling regime. By in- 
creasing the imbalance progressively, we expect the emer- 
gence of angular dependent magnetization distribution, 
indicating the existence of angular FFLO state. Very re- 
cently, the MIT group produced preliminary experimen- 
tal evidence for superfluidity of ultracold ^Li atoms in 
optical lattice We expect that a future experiment 
for imbalanced fermions will be conducted to examine 
our prediction. An alternative scenario to detect such 
an exotic state is to carry out the experiment in a thin 
superconducting ring of heavy fermion materials. By ap- 
plying strong magnetic field parallel to the plane, the 
angular-dependent FFLO state may show up due to the 
special topology of the ring structure. Probe of the real 
space modulation of the pairing gap as well as magne- 
tization can be achieved by using SQUID or STM tech- 
niques. The third proposal is to generate a Mexican hat 
trapping potential in 3D so that the distribution of the 
confined fermionic atoms may form a donut-like struc- 
ture. In such a case, the angular dependent distribution 
of magnetization profile can be directly measured by us- 
ing the time-of-flight imaging techniques. 

In summary, we have explored theoretically the novel 
superfluidity and ferromagnetism in ultracold fermionic 
atoms on a 2D lattice combined with a harmonic trap. 
Due to the interplay between the confined ferromagnetic 
and superfiuid orders, exotic FFLO phases have been re- 
vealed. At low density, our theory shows a bimodal distri- 
bution of the magnetization profile and a fermion pairing 
gap oscillating along the radial direction, reproducing the 
main features observed in recent experiments. We pre- 
dict a more exotic angular FFLO state at high densities 
and square lattice like FFLO state at the intermediate 
densities. 
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